"Unusual" critical states in type-II superconductors 
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We give a theoretical description of the general critical states in which the critical currents in type- 
II superconductors are not perpendicular to the local magnetic induction. Such states frequently 
occur in real situations, e.g., when the sample shape is not sufficiently symmetric or the direction 
of the external magnetic field changes in some complex way. Our study is restricted to the states in 
which flux-line cutting does not occur. The properties of such general critical states can essentially 
differ from the well-known properties of the usual Bean critical states. To illustrate our approach, 
we analyze several examples. In particular, we consider the critical states in a slab placed in a 
uniform perpendicular magnetic field and to which two components of the in-plane magnetic field 
are then applied successively. We also analyze the critical states in a long thin strip placed in a 
perpendicular magnetic field which then is tilted towards the axis of the strip. 

PACS numbers: 74.25.Sv, 74.25.Qt 



I. INTRODUCTION 



The concept of the critical state introduced by Charles 
Beani. is widely used to describe various physical phe- 
nomena in the vortex phase of type-II superconductors, 
see, e.g., Refs. [HH and citations therein. According to 
Bean, in the critical state of type-II superconductors with 
flux-line pinning, the driving force of the currents flowing 
in this state is balanced by the pinning force acting on 
the vortices. The critical state is characterized by the 
component of the current density flowing perpendicular 
to the flux lines, j c ±, since only this component generates 
a driving force. It is assumed in the critical-state theory 
that this j c ± is known, i.e., it is a given function of the 
magnetic induction B, j c j_ = j c ±(B), and the problem of 
this theory is to find the appropriate distribution of the 
magnetic fields and currents in the critical state. Below, 
for simplicity, we shall assume that the magnetic fields 
H in the superconductor considerably exceed the lower 
critical field H c \, and so we put B = /j,qH throughout 
the paper. Beside this, we deal only with bulk super- 
conducting samples, assuming that all their dimensions 
noticeably exceed the London penetration depth, and we 
consider the critical state macroscopically, averaging vor- 
tex structures and the appropriate microscopic currents 
over a scale exceeding the intervortex spacing. 

Hereafter we shall call the critical states "Bean critical 
states" if the current density j is perpendicular to the 
local magnetic field H at every point of a superconductor, 
j = j_L) and thus j = j± = j c ±. This definition imposes 
limitations on the direction of the currents in the critical 
state, but it does not imply constancy of j c ±, e.g., j c ±{H) 
can be as in the Kim models The Bean critical states can 
be found from the static Maxwell equations, 



rotH=j, divH = 0, 



(1) 



and the conditions on the current density 

divj = 0, (2) 

3± = 3c±, j|| = 0, (3) 

where ju is the component of the current density along 
the local magnetic field H. Such states usually occur 
when the shape of the superconductor is sufficiently sym- 
metric and the external magnetic field H a is applied 
along a symmetry axis, so that the direction of the cur- 
rents is dictated by the symmetry of the problem. Most 
of the known solutions of the critical state problem de- 
scribe just these Bean states. For example, this is the 
well-known solution for an infinite slab in an external 
magnetic field parallel to its surface^ and also the solu- 
tion for an infinitely long cylinder with arbitrary cross- 
section in a magnetic field parallel to its axis since the 
currents flow perpendicular to this axis* 2 Bean critical 
states also occur in infinitely long and thin strips^&i 
and in thin disks- in a perpendicular magnetic field even 
if j c ± depends on |B| = B or on the angle between B 
and the normal to the sample plane d 10 ' 11 If the applied 
magnetic field is tilted to the plane of an infinitely long 
strip 12 ' 13 i 14 i 15 or slahlS but remains perpendicular to the 
sample axis, the critical currents flow along this axis, 
and a Bean critical state occurs. Further examples of the 
Bean critical states in samples of a complex shape can 
be found in Ref. U3IIIGJ[ioHl A characteristic feature 
of all these Bean critical states is that the perturbation 
of the current distribution caused by a change of the ap- 
plied field propagates into the sample as a sharp front at 
which the direction of the currents changes abruptly. 

In real samples of nonsymmetric shape, or when the 
applied magnetic field changes not only in amplitude but 
also in its direction, adjacent flux lines may be slightly 
rotated relative to each other in the critical state. This 
rotation generates a component of the current along the 
magnetic field^ jii . The rotation of flux-lines can lead to 
their mutual cutting! 2 ^ Flux line cutting occurs when 
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the component of the current density parallel to the mag- 
netic field, ju , exceeds some longitudinal critical current 
density j c u. In this situation a vorte?*^ or a vortex 
array 2 ^ becomes unstable with respect to a helical distor- 
tion, and the growth of this distortion leads to flux-line 
cutting. When both ju and j c ± are equal to their crit- 
ical values j c u and j c ±, respectively, the so-called dou- 
ble critical stat o 22 ' 25 occurs in the superconductor.— For 
example, this state appears in some region of a super- 
conducting sampl o 25 ! 27 when a rotating magnetic field of 
constant magnitude is applied to a superconducting disk 
(or slab) in its plane. 28 ' 29 ' 30 The double critical state can 
be still described by Eq. |T|), but with the following 
conditions on the current density j = jj_ + jy: 

jx = jcx, j|| = 3c\\- (4) 



As was pointed out in Refs. [25ll27l . one more type of 
critical states can exist in superconductors. In these 
states j± < j c j_ and ju = j c u, i.e., only flux cutting oc- 
curs without any transport of vortices. The description 
of such C-critical states (C means cutting) in samples of 
arbitrary shape can be obtained by an immediate gener- 
alization of the approach used in Ref. [25| 27 for a super- 
conducting slab. Below we shall not analyze such states 
in detail but only briefly outline this generalization. 

In Sec. II of this paper we develop the approach of 
Ref. 1371 - In particular, we take into account the de- 
pendence of j c ± on ju and anisotropy of flux-line pin- 
ning. We also discuss the relationship between the equa- 
tions of Ref. [37] and the variational principle recently 
proposed i 44 ! 45 ' 46 In Sec. Ill we then analyze three ex- 
amples of the general T-critical state. 



The concept of the critical state with flux-line cutting 
was further developed in Refs. 13^13311341 to explain the 
observed suppression of the magnetic moment of a su- 
perconducting slab under the action of an ac magnetic 
field^i^ 6 . 

However, in many real situations a change of the di- 
rection of the external magnetic field or a nonsymmetric 
shape of the sample does not lead to flux cutting in the 
superconductor, i.e., ju does not reach j c n in the critical 
state. In such situations there is no explicit condition on 
the magnitude of j|| except that ju < j c n, and the static 
equations (P) and ^ with the only restriction j± = j c ± 
are not sufficient to find the distributions of the mag- 
netic field H(r) and current density j(r) in the critical 
state. This problem for the special case of a slab with 
an in-plane magnetic field was solved in Refs. 25,27. The 
full set of the critical-state equations for arbitrary shape 
of the sample and for any quasistatic evolution of the 
vector of the applied magnetic field H a was obtained in 
Ref. H3, where it was also shown that in contrast to the 
common Bean critical states, a perturbation of the cur- 
rent distribution in such critical states propagates into 
the sample smoothly in a diffusive way. We emphasize 
that this class of critical states with ju < j c || corresponds 
to the general situation, while the common Bean criti- 
cal states and the double critical states are only limiting 
cases occurring when jy = or jy = j c ^ , respectively. 

Such general critical states, which we shall call the T- 
critical states (T means transport) j 3 - 8 - occur even for sim- 
ple experimental situations. In particular, they appear in 
a certain region of thin rectangular platelets in a perpen- 
dicular magnetic field (in platelets with thickness exceed- 
ing the London penetration depth this is the region which 
is not penetrated by the perpendicular component of the 
magnetic field) . 3 - Critical states of this type also appear 
at the vortex-shaking in rectangular platelets 4 ^ and even 
in strips if the ac field is along the axis of the strips^ 
They also occur in low-frequency ac experiments with a 
slab when a circularly polarized ac field is applied per- 
pendicularly to the dc magnetic field H a that is normal 
to the plane of the slab i 42 ' 43 



II. GENERAL CRITICAL STATES 

A. Critical-state equations 

The critical state is well established in a sample if the 
characteristic time of change of the applied magnetic field 
H a , j c j_d/\dH a /dt\, considerably exceeds the time of flux 
flow across the sample, pgd 2 /ps, where d is a character- 
istic size of the sample and pg is the flux-flow resistivity. 
In other words, the concept of the critical state can be 
used for a description of the magnetic-field and current 
distributions in superconductors if the generated eddy 
electric fields are relatively small, 



Hod 



dt 



< PsjcX- 



(5) 



The ideal critical state thus corresponds to the limit 
ps — > oo. Below we imply condition © to be fulfilled. 

The general T-critical states with ju < j c || can be de- 
scribed by the following approach: 37 The static equations 
(fTJ) and ^ are supplemented by the quasistatic Maxwell 
equation 



rot E = — /ioH, 



(6) 



where H = d~H/dt, and E is the electric field gener- 
ated by a change of the applied field H a . For the set of 
equations (fTJ), ([2]), and ([6]) to be solvable, it has to be 
supplemented by the current-voltage law E^B)^ This 
law is introduced from two well-known physical ideas: 1) 
At any given j and B, the direction of E follows from 
E = [B x v], i.e., 



E|| [Bxv], 



(7) 



where v is the vortex velocity caused by the Lorentz force 
[j x B]. Here for simplicity we shall neglect the so-called 
Hall angle^ and so the directions of v and the Lorentz 
force coincide. 2) The magnitude of E is found from the 
condition that 



LLl 



: ]cX- 



(8) 
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In fact, this condition may be interpreted as the following 
current-voltage dependence: 



|E| = at j x < j c x, 
|E| —> oo at j± > j c ±, 



(9) 



which just corresponds to the ideal critical state. 

To proceed with our analysis, let us introduce the fol- 
lowing notations for the magnetic field H(r) and the cur- 
rent density j(r) in the critical state: H(r) = H(r)v(r), 
j(r) = j(r)n(r) where H and j are the absolute values 
of the magnetic field and the current density while the 
unit vectors v and n define their directions. Then, the 
component of the current density perpendicular to the 
magnetic field is given by 

3± = j - ^.j) = jc±n±(r). 

Here the unit vector ni defines the direction of 
ni = (n — i/(i/n))/ D; D = \/T— (nv>) 2 is the normaliz- 
ing factor that is equal to the sine of the angle between 
H and j, and we have taken into account the condition 
= j c ±. These formulas also lead to the explicit ex- 
pression for the magnitude j of the current density, 



J = 



Jc ±(H) 
D 



(10) 



that is only another form of the condition = j c ±- Let 
us now formulate condition ([7]), Let at a moment of time 
t the external magnetic field H a (t) change infinitcsimally 
by 5H a = il a 8t. Under the change of H a , the critical 
currents locally shift the vortices in the direction of the 
Lorentz force [j x v]; this shift generates an electric field 
directed along \v x [j x u}] = i.e., along the vector 
ni. Thus, we can represent the electric field E(r) in the 
form: 



E 



(11) 



where the scalar function e(r) is the modulus of the elec- 
tric field. Note that the electric field in general is not 
parallel to the total current density j (r) . With formulas 
(fTO)) and CLI]), equations ©, ©, and ([6]) are sufficient 
to describe the T-critical states in a sample of arbitrary 
shape. It is important that the magnetic fields H(r) and 
currents j(r) in the critical state at the moment of time 
t + St depends only on the field and current distributions 
in the previous critical state at the moment t and on the 
change of the external field <5H a = H a <5<, while the elec- 
tric field e is proportional to the sweep rate H a rather 
than to 5H a , and so it plays an auxiliary role in solving 
the critical-state problem^ 

We emphasize that e is now found as a solution of the 
set of equations H}, J2|), ©, (10]), (LI} without using 
the specific current-voltage dependence ©. The explicit 
equation for the scalar function e(r) has the formal 



n_L- {rotrot(E) - (v ■ rotv) rot(E)} 
djcx(H) 



<9H 



rot(E), (12) 



where E is given by Eq. (|11[) . Continuity of the magnetic 
field on the surface of the superconductor, S, yields the 
boundary condition to Eq. (fT2"|) : 



rot (E(r s )) 



Mo 



[R x rotrot(E(r'))] 
4ttR 3 



dr',(13) 



where r$ is a point on the surface S, R = rg — r', 
R = |R|, and the integration is carried out over the vol- 
ume of the sample. The right hand side of this boundary 
condition expresses fioil on the surface of the supercon- 
ductor (but reaching from outside) with the use of the 
Biot-Savart law. If in the critical state of the supercon- 
ductor there are also boundaries at which the direction 
of the critical currents changes discontinuously or which 
separate regions with j± = j c ± from regions with j = 0^ 
the function e(r) has to vanish at these boundaries to 
provide continuity of the electric field enj_ there. 

In practical calculations of critical states developing in 
the process of changing H a (t) it is convenient to rewrite 
Eqs. |lj and ((6]) in the form 



^oj[ n x h] = — [n x rotrot(E) 



(14) 



which is a differential equation for the angles defining 
the direction of j, i.e., the unit vector n = j/j. Note that 
since the distributions of the magnetic fields and currents 
in the critical states of a superconductor are independent 
of the sweep rate H a , their temporal dependence is only 
a parameterization of their dependence on H a . 

Let us now write explicitly the applicability condition 
of the above theory. Since the projection of j on the local 
direction of H is j c ±(nv) / D, the condition that flux- line 
cutting is absent leads to the following restriction on the 
angle between the local j and H: 



v/i - (^) ! 



< 



Jc|| 
jc±' 



(15) 



where j c \\ is the longitudinal critical current density. 

Finally, we make several remarks on the electric field. 
It may turn out that the electric field en± obtained with 
Eq. (fP2)) does not satisfy the condition div(enjj = 0. To 
clarify this situation, it is necessary to remember that a 
moving vortex generates an electric dipolar moment, 48 
and hence the moving vortex medium is characterized 
by the vector of polarization P which is the macroscopic 
density of this moment. It follows from the results of 
Ref. [48| that P = -eni, and a nonzero div(enjj means 
that in a type-II superconductor the electric-charge den- 
sity — divP appears which generates a curl-free electric 
field E p = — V$ described by the scalar potential 
This potential field is a part of the total electric field 
given by E = enj_ inside the sample, and it obeys the 
equation divE p = div(en^), i.e., 



A$ = — div(enx), 



(16) 



where A = divV. At the surface of the sample, S, the 
field E p satisfies the same boundary conditions as in the 
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electrostatics of dielectrics 4i The tangential components 
of E p and the normal component of E p + P = E p — en± 
are continuous there. Since P = outside the sample, 
the latter condition means that 

(E+-E;)T = -en ± T (17) 

where E+ and E~ are the surface potential fields calcu- 
lated outside and inside S, respectively, and r is the nor- 
mal to S pointing outside. The right hand side of Eq. (fT7|) 
gives the surface-charge density induced by moving vor- 
tices in the sample. Note that the potential part of eni 
does not influence the magnetic fields and currents in the 
critical state since rotE p = 0. Appearance of this part 
is caused by condition (fTTj) that dictates the direction of 
the electric field in the sample. Although both the induc- 
tive part of the electric field, en± — E p , which generates 
the critical states, and the potential part E p can be mea- 
sured in certain situations^ we shall not analyze electric 
fields in detail in this paper since these fields plays only 
an auxiliary role in the critical state problem. See also 
the recent book on electric fields £i 

Generally speaking, in the process of changing H Q 
a migration of the induced charges p = div(en^) oc- 
curs, which leads to a generation of currents satisfying 
divj = —(dp/dt) and violating Eq. (J2J). However, these 
nonstationary currents are proportional to the second 
power of the sweep rate H a and are negligible under as- 
sumption ((5|). 

B. Generalizations 

We now point out some generalizations of the above 
results which may be useful in analyzing critical states in 
real situations. 



1. j c ± depends on j\\ 

The current-voltage law used in Sec. II A, Eq. ([5]), 
means that flux creep is negligible in our approach. In 
this case the critical current density is found from the 
condition that the creep activation barrier U of a vor- 
tex bundle is equal to zero. It has been implied above 
that j c _i_ may depend on B but is completely indepen- 
dent of the magnitude of j\\ . In other words, the form 
U = U(j±,H) has been assumed for this U. How- 
ever, the creep activation barrier U, generally speak- 
ing, may depend not only on j± and B but also on the 
j'll that characterizes flux-line misalignment in the bun- 
dle, i.e., in the general case one has U — U(j±, j||,B). 
Then the critical current density j c ±_ determined from 
U(j ±i J||,B) = takes the form j c± = j c± (B,j\\). One 
may expect that this dependence of j c ± on the longitudi- 
nal current component j\\ is especially noticeable when ju 
is close to its critical value j c », and hence j c ±(B, j c \\) m 
general differs from j c ±(B,0). Similarly, the activation 



barrier U cn t for flux cutting is a function of both current- 
density components and of the magnetic induction, i.e., 
Ucut = U cut (j±,j\\,B), and the condition U cut = gives 
3c\\ = Jell (B, j±). In Fig. la at a fixed B we schematically 
show the dependences of j c ± on ju and of j c u on j± in the 
plane (i_i_,i||). Note that these dependences cross when 
the equations U(j±,j»,B) = and U cut (j±, j\\ , B) = 
hold simultaneously. This occurs at isolated points in 
the plane since the barriers U and [7 cut character- 

ize different physical processes and are essentially differ- 
ent functions of the current components. These points 
correspond to the double critical states when j± = j c ± 
and j|| = j c n. In Fig. la the top/bottom and right/left 
sections of the curves between the four points describe 
jc±(j\\) fii the general T-critical state and j c \\{j±) i n the 
C-critical state. 

The dependence j c ±(j\\) leads to a replacement of 
j c ±(H) by j' c _l(H, j\\) in formula (I10| that now reads 

jD = j c± (H,jy/l-D 2 ). (18) 

The dependence j c ±(j\\) also leads to a modification of 
Eq. (fT2| . In the right hand side of this equation the term 
—fj,a(dj c ±/dj\\)(dj c \\/dt) should be added that equals 

%^ ( % n x ■ rot(E) + v ■ rot rot(E) ) (19) 

with E from Eq. pT|) . Note that the first term in this 
expression has no singularity at H — > since the combi- 
nation n^- rot(en^) can be also rewritten as enj_- rotn^ 
and e = | [B x v] | cx H. 

In Refs. I52II53I a phenomenological model of the general 
critical state was considered that described sufficiently 
well a number of experimental data on the magnetization 
of a slab and of a disk in magnetic fields parallel to their 
planes. In fact, in this model a certain type of the de- 
pendence of j c ± on ju (and of j c ii on j±) was introduced. 
Even though the directions of the electric field in this 
model do not satisfy the physical requirement (jTTJ) , the 
sufficiently good description of the data seems to indicate 
the importance of this dependence in real situations. 

2. Anisotropy of flux-line pinning 

In deriving Eq. (|11[) we have assumed that when H Q 
changes, vortices shift in the direction of the local Lorentz 
force [j x B]. However, in the case of anisotropic pinning 
this assumption may fail. Nevertheless, even in this case 
the direction of the shift can be expressed via the direc- 
tions of j and v = H/iJ, see Appendix A in Ref.[39|. Now 
the unit vector u along the electric field, E = ue, is 

u = nj_ cos 6 + [u x rijj sin S, (20) 

where the angle 8 describes the change of the direction 
of the electric field due to anisotropic pinning. If in the 
plain perpendicular to the local H the critical current 
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current-voltage dependence: 
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FIG. 1: (a): Dependences of j c ± on jy and of j c y on j± (solid 
lines) shown schematically in the j±-j\\ plane at a fixed B, 
see Sec. II Bl. The crossing points of the lines correspond to 
the double critical states when j± — j c ± and ju = j c ii. Shown 
are also the directions of the electric field for the appropriate 
critical states, (b) and (c): The boundaries F (solid lines) 
of the regions A introduced by Badi'a and Lope a 44 ' 45 ' 46 in 
the j±-j\\ plane, see Sec. II C. Here the regions A are (b) a 
circle and (c) a rectangle. Shown are also the directions of the 
electric field according to the approach of Badi'a and Lopez. 



density j c ± depends on its direction n^, the angle 5 is 
found from^S 



tan 5 = — 



0[lnj e± (B,, 

dcj> 



(21) 



where <fi is the angle defining the direction of = 
(cos </>, sin (f>) in the plane perpendicular to H. When j c ± 
is isotropic in this plane, we obtain 5 — 0, and thus u 
coincides with n^. 

Equations (|20|) and (|2Tjl give a relation between nj_ 
and u. When S ^ 0, i.e., when the vector u differs from 
, the only change in the critical state equations is that 
eni in Eq. (fTTj) is replaced by eu, and j c ± (H) in Eq. (flQ|) 
is now i c ±(H, n ± ). 



3. C-critical states 

As it was mentioned in the Introduction, in the case 
of an infinite slab the critical states with flux-line cut- 
ting but without flux-line transport were considered in 
Refs. [25ll27l . For samples of arbitrary shape such C- 
critical states can be described by Eqs. {U, ©> and 
(|6|), but now the electric field is along the local H, i.e., 
E = ve. This condition replaces Eq. [or (|TT))]. The 
absolute value e of the electric field is now determined by 



|E| = 
|E| -> 00 

and leads to the formula 



at in < icih 

at i|| > i c ||, 



icii(H) 



nu 



(22) 



(23) 



Equations 
tively. 



replace Eqs. §§§ and (JTUJ), respec- 



C. Variational principle 

Recently ; 44 ' 45 i 46 a variational principle was put forward 
to describe the critical states in superconductors. In de- 
riving this principle Badi'a and Lopez used Eqs. (Q]), © 
and the current-voltage law with |E| = when j is inside 
some region A of the j-space and |E| — * 00 when j lies 
outside this region. In other words, the critical states cor- 
respond to the boundary T of the region A, see Fig. lb,c. 
However, the physical idea of the direction of the electric 
field, Eq. ([7]), was not incorporated in their principle. In- 
stead of this they find the direction of the electric field 
from some condition of maximality of their Hamiltonian. 
This leads them to the conclusion that the electric fields 
in the critical states are directed along the normals to 
the boundary T at the appropriate points, Fig. lb,c. 

Within our approach their boundary V corresponds to 
the contour composed of the dependences i c j_(i||) and 
ic||(ij-); see Fig- la and Sec. II Bl. But in our general 
T-critical states with jn < j c \\ the electric field is always 
perpendicular to the local H (i.e., to jy), and in the C- 
critical states with flux-line cutting but without flux-line 
transport the electric field is along the local H. It is 
clear that only in the case when A is a rectangle does the 
approach of Badia and Lopez lead to the correct results 
for the electric field, Fig. lM- However, in general their 
approach leads to contradiction with existing physical 
concepts ^ 23 ' 24 In particular, in the so-called isotropic 
model, when A is a circle, Fig. lb, the electric field E 
is parallel to j, and hence a nonzero E along H appears 
even for an infinitesimally small longitudinal component 
of j, i.e., flux- line cutting in that model occurs without 
any threshold j c \\ ^ 



III. EXAMPLES 

We first consider two examples of the general critical 
state in an infinite slab of thickness d. Let this slab fill 
the space |x|, \y\ < 00, \z\ < d/2, and be in a constant 
and uniform external magnetic field H az directed along 
the z axis, i.e., perpendicularly to the slab plane. The 
critical current density j c ± is assumed to be constant 
in this slab. In the first example a constant field H ax 



G 



(H ax > J c /2 = dj c ±/2) is applied along the x axis, and 
after that the magnetic field H ay is switched on in the y 
direction. This example was considered in our paperj2Z 
but there H ax , H ay , and J c were assumed to be small as 
compared with H az , i.e., the tilt angle 9 of the magnetic 
field to the z axis was always small. Now we do not 
put this restriction, and the angle 9 may be sufficiently 
large. But we still assume that flux-line cutting does not 
occur (see below). This example may be considered as 
a modification of the experimental conditions of Refs. 
where the suppression of the magnetic moment of the 
slab was investigated at H az = 0. In the second example 
the critical current along the y axis is applied to a slab, 
and after that the magnetic field H ay is switched on in 
the same direction. 

The critical state equations are the same for these two 
examples. The difference is only in the boundary con- 
ditions. Let us write these equations. The condition 
divj = dj z /dz = together with j z \ z =±d/2 = yields 
j z = 0, i.e., the currents flow in the x-y planes^ Then, 
to describe the critical state, we may use the parameter- 
ization: 



n± y 
n± z 



sin ip — sin 2 9 sin ip cos(<p — ip) 
D : 
sin 9 cos 9 cos(p — ip) 
D ' 



(28) 



and equation (fT2")l for the electric field e takes the form: 

n± x (en± x )"+ n ±y (en ±y )" 

- ip 1 sin 2 9(n' ±x n± y - n ±x n' ±y )e = 0, (29) 

where the dash over a symbol means d/dz. For the angle 
ip we obtain from Eq. (Q2 



jc± dp . , .// . /on\ 

^ a ~D~dt = v cosip - (enj_ x ) smp. (30) 

At small 9 when D 1 and « n, equations (12T)|) and 
(|3H)) reduce to the form that was used in Ref. [33 



"-{p'fe = 0, 

dip . . ,, 

MoJ c ±^- = 2e ip + eip . 
at 



(31) 
(32) 



j =j c ((p,9,ip)(cosip(z), sin<^(z),0), 
H(z) = zH az + h(z), 
h(z) = (h x (z),h y (z),0), 



(24) 



where z is the unit vector along the z axis; j c ((p,9,ij)) 
is the magnitude of the critical current density when 
a flux-line element is given by the angles ip and 9, 
tamp — h y /h x , tan# = (/i 2 + hy) 1 ^ 2 / 'H az , while the cur- 
rent flows in the direction defined by the angle ip; all these 
angles generally depend on z. A dependence of j c on the 
orientation of the local H appears even at a constant j c ± 
ii j c is not perpendicular to this H, and this dependence 
is described by formula (fTU]) , where D in terms of the 
angles is 



D=[l- cos 2 (ip - ip) sin 2 9] 1/2 . 



(25) 



Equations (f2T>)) - (|30[) provide the complete description 
of any general T-critical state in an infinite slab with 
j c j_(H) = constant in the absence of flux-line cutting. 
For different critical-state problems only the boundary 
conditions should be appropriately chosen. Note that 
the usual Bean critical states in the slab correspond to 
discontinuous solutions tp(z) of these equations. 

In the case of the slab condition (fT5"f of absence of 
flux-line cutting leads to the following restriction on the 
angles 6, ip, and ip: 



sin0| cos(<£ — ip)\ j c \\ 



(33) 



This condition is fulfilled at any ip and ip, i.e., at any 
direction of j and h, if the z component of the magnetic 
field, H az , is not too small, 



With this parameterization, the equation divH = is 
satisfied identically, while the Maxwell equation rotH = j 
reads 



dh x jc± sin ip 
~dz~ _ D ' 
dhy jci. cos ip 



dz 



D 



(26) 
(27) 



These equations differ from the appropriate equations of 
Ref. [37] by the factor 1/D, which is not unity now. 

In the case under study one has n = (cos <p, sin ip, 0), 
v = (sin 9 cos ip, sin 9 sin ip, cos 9). Then, a direct calcu- 
lation gives the following expressions for the vector nj_ 
defining the direction of the current component perpen- 
dicular to H: 



cos <p — sin 2 9 cos ip cos(<p — ip) 
D ! 



I hi 4. hi 
tan 9 = — < — - 



H„ 



Jc± 



(34) 



We imply this condition to be fulfilled below. 



A. First example: H ax and H ay 

In the first example that we consider, a constant field 
H ax (H ax > J c /2 — dj c ±/2) is applied along the x axis, 
and after that the magnetic field H ay is switched on in the 
y direction. Then, the boundary conditions to Eqs. (f2"6"|) 
- ((30]) at z = d/2 are 



(en± x )' = -/i 



h x — H ax , h y — H ay (t), (35) 

dH av (t) 



dt 



(en± v )' = 0, 



(36) 
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or equivalently, conditions (|36|) . which follow from for- 
mula (|13p . can be rewritten in the form: 



/ / / / dH a y (i) f 
e [n ±x nx y - n± x n ±y ) = Mo n 



e(n' ±x nj_ y - n± x n' x ) = -/i 



dt 

dH ay (t) 
dt 



(37) 



Taking into account the symmetry of the problem, 
e(— z) = e(z), ip(— z) = ip(z) — n, h(— z) = h(z), it 
is sufficient to solve equations (|2l))) -(l3T)l) in the region 
< z < d/2. At z = 0, where the direction of the 
currents changes discontinuously, one has the additional 
condition for e, 



e(0) = 0. 



(38) 



Since after switching on H ax , the critical currents flow 
in the y direction, we have the following initial condition 
for Eq. (|30]): 



ip(z,t = 0)=7r/2, 



(39) 



where the moment t = corresponds to the beginning of 
switching on H ay . As to the initial magnetic-field pro- 
files, equations ([27]), and (32]) give h y (z,t = 0) = 
H a y — and h x {z, t = 0) = H ax — 0.5 J c + j c ±z where 
J c = jc±d. 

In the limiting case H az 3> H ax , H ay , J c , the solution 
of equations (j2l)]) - (f3T)l) with conditions ([3"5]) - (|3^|) was in- 
vestigated in Ref. l37|. Since in this case w n, one 
finds that the electric field en± is along the current den- 
sity jn, and in fact, we arrive at a situation which can 
be formally described by the so-called isotropic model of 
Badfa and L6pez4£ As was explained in Sec. II C, this 
model in general does not lead to the correct direction of 
the electric field. In particular, it fails in the following 
situation discussed by Badfa and Lopez 4^ A slab with 
H az — 0, H ax — constant, and oscillating H ay . But in 
the case H az ^> H aX} H ay , J c , which in reality was not 



considered in Ref. |45j, the isotropic model gives the cor- 
rect results, and the numerical data of Ref. |45| agreed 
with those of Ref. [37| and can be used to describe this 
limiting situation. 

In the case H az ~ H ax , H ay , J c , the solution of equa- 
tions ([25])- 1301) with conditions (|55|) - (f3"9"]) is shown in 
Fig. 2. We present (p(z), ip(z), h x (z), h y {z) in the se- 
quence of the critical states developed in the process of 
increasing H ay . We do not show the electric field e that is 
proportional to the sweep rate H ay and plays an auxiliary 
role. As was noticed previously^ in stark contrast to the 
Bean critical states, in which any change of the current 
direction occurs inside a narrow front, in the general T- 
critical state the change of the angle f(z) with increasing 
H ay has diffusive character. But there is a difference be- 
tween the data of Fig. 2 and the results^, obtained in the 
case H az 3> H ax , H ay , J c when the currents in the criti- 
cal state are almost perpendicular to the local magnetic 
fields. In the latter case at H ay > J c the field profile 



^ 0.6 
0.4 




_2_ . 
5.4. 



H =4 H =1.1, H =1.5, diamag 




FIG. 2: Top: Profiles of the angle of the currents ip(z) 
(solid lines), and of the field angle il}(z), tan^ = h y (z)/h x (z) 
(dashed lines) in the critical states of the slab at H a z = 1.5, 
H ax = 1.1 and H ay = 0.01, 0.05, 0.1, 0.2, 0.35, 0.5, 0.7, 1, 
1.4, 2, 3, 4, 6, 9. Bottom: The magnetic field components 
h x (z) (dashed lines) and h y (z) (solid lines) in the same criti- 
cal states. We start from the diamagnetic initial critical state 
with h x (z) = H ax — 1 + 2!, h y (z) = H ay = 0, and ip(z) = 7r/2. 
Here z is in units of d/2, and the magnetic fields in units of 
j c± d/2 = J c /2. 



h x (z) becomes practically constant and coincides with 
H ax , while the angle <p tends to 7r. On the other hand, 
we sec from Fig. 2 that for H az ~ J c the angle ip lies in the 
interval tt < ip < 3ir/2 at H ay > J c . In other words, the 
y component of j(z) has the opposite direction as com- 
pared with the initial state. This leads to the fact that 
at H ay > J c the field h x increases towards the central 

plane of the slab, z = (but h(z) = ^Jh x + h y still de- 
creases towards this plane), and the initial diamagnetic 
state with the magnetic moment M x = —j c ±d 2 /4 (per 
unit area) turns into a paramagnetic state with positive 
M x . 

In Fig. 3 we show the same sequence of the critical 
states but in the case of the paramagnetic initial state. 
This initial state is obtained if one first increases the field 
H ax essentially above the field of full flux penetration and 
then decreases it to a prescribed value. Now the initial 







0.2 



0.4 0.6 



1 2 3 , , 4 5 6 7 



FIG. 3: The same as Fig. 2, but at H ay — we start from the 
paramagnetic initial critical state with h x (z) = H ax + 1 — z, 
hy(z) — 0, and f(z) = — 7r/2. For clarity the profiles <p(z) 
and h x (z) at Zf aa > 2 are shown by dash-dotted lines (and 
for H ay < 1.4 by solid lines). 



condition to Eq. ((501) i s 

i = 0) = -tt/2, (40) 

and the magnetic fields at t = are given by h y (z,t — 
0) = H ay = 0, M-M = 0) = H ax +0.5J c -j cA _z. It is seen 
from Fig. 3 that although at H ay < H az a decay of the 
initial paramagnetic profile h x (z) occurs, with a further 
increase of H ay new paramagnetic states are developed 
that are close to the appropriate states of Fig. 2. 

In Fig. 4 we compare the if^-dependences of the mag- 
netic moment (M x , M y ) per unit area of the slab, 58 

fd/2 

M = / z[z xj]dz, (41) 

obtained using the two sequences of the profiles j(z, H ay ) 
developed from the diamagnetic and paramagnetic ini- 
tial states with the same H ax . It is seen that in both 

cases \M X \ and \M y \, and even more M = + M£, 

can exceed the "saturation value" j c ±d 2 /4 used as unit in 
Fig. 4. This is possible since the current density j exceeds 
j c ± when it does not flow at a right angle to the vortices. 



FIG. 4: The magnetic moment (M x ,M y ) of the slab per unit 
area denned by Eq. (|4ip as a function of H ay . Shown are the 
diamagnetic (solid lines, see Fig. 2) and paramagnetic (dashed 
lines, see Fig. 3) cases for H ax = 1.1 and H az ~ 1.5 (top) and 
Haz = 10 (bottom). Units are j c ±d/2 for H and j c ±d 2 /4 for 
M. 



This excess of j leads to that M x (H ay ) does not saturate 
at large H ay but continues to increase nearly linearly, 
with slightly negative curvature. The other component, 
M y (H ay ), at large H ay practically saturates to a value 
slightly lower than — j c j_d 2 /4. Of course, one should keep 
in mind that in reality the region of large H ay where these 
results for M are applicable is limited by condition f|34[) . 
Note also that in agreement with Figs. 2 and 3 the mag- 
netic moment M x (H ay ) is always positive at sufficiently 
large H ay , and the diamagnetic and paramagnetic initial 
states lead to practically the same M x at such H ay . 

As it is known, field-cooled type-II superconducting 
samples frequently exhibit a positive magnetic moment; 
see, e.g., paper— and references therein. Different expla- 
nations of this paramagnetic effect were put forward. In 
particular, this effect may be associated with the com- 
pression of trapped magnetic flux in the sampled The 
data of Figs. 2-4 show that in principle, the paramagnetic 
effect may be also due to the field-cooling caused gener- 
ation of critical states in which the circulating currents 
are not perpendicular to the local magnetic fields. 

The general T-critical states considered here can be re- 
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FIG. 5: Profiles of the angles (/'(z) (solid lines) and tp{z) 
(dashed lines) (top), and of the magnetic field components 
h x (z) (dashed lines) and h y (z) (solid lines) (bottom) for a slab 
with applied sheet current J = J c /2 and H ax = 0. The pro- 
files are calculated from Eqs. (|43)) and (gHJ for H ay = 0.796, 
0.81, 0.85, 1, 1.4, 2, 3, 5, 10, 20, 40, 100, 500 and large 
H az = 20 (all the magnetic fields are in unit of J c /2). 
Hay = 0.796 is the field of full flux penetration at J = J c /2. 



B. Second example: J v and H ay 

We now consider the second example of the general T- 
critical state in the slab. It is assumed that the slab is in 
a uniform magnetic field H az along the z axis, the current 
J (per unit length along x) flows in the y direction, and 
at t — the field H ay is switched on. In this case the 
boundary conditions at z — d/2 are 



J 



(enxx)' = -Mo 



dH ay (t) 
dt ' 



-^ay (^) ) 

{en± y Y = 0. 



(42) 
(43) 



The symmetry of the problem is now described by the 
relationships: e(—z) = e(z), tp(—z) = Tr-ip(z), h x (—z) = 
—h x {z), h y (—z) = hy[z), and at z = the direction of 
the critical currents changes continuously. Thus, instead 
of condition (|3"51) we have at z = 0, 



M0) = 0, <p(0)=tt/2. 



(44) 



As in the first example, we shall consider the critical 
states only in the interval < z < d/2. 

If J is less than J c — j c ±d, in the initial state the 
current flows only at 1 — (J/J c ) < 2z/d < 1. After 
switching on H ay the current distribution develops over 
the whole thickness d when H ay reaches a penetration 
field H® y < J c /2, and we shall analyze the critical states 
only after this penetration of the current has occurred, 
i.e., at H ay > H° y . 

Below we consider only the case H az ^ Jc- In this 
case in the leading order in the small parameter J c /H az 
we find the following analytic solution of Eqs. (|26|) - (f3"U|) 
with boundary conditions (|42|l - (l44|) . see also Fig. 5: 



Ho H ay a D 
sin ip cos 9 ' 

cos 8 \ a 



cot ip = — - 



z cos 9 



(45) 



alized in experiments similar to the experiments of Park 
et al&L and Fisher et al.^ except that now the field H az 
perpendicular to the plane of the sample is not equal to 
zero. Such investigations would enable one to compare 
the theoretical results for the general T-critical states 
with the appropriate experimental data avoiding compli- 
cations due to flux-line cutting. To prepare the initial 
state which was described above, e.g., in a supercon- 
ducting strip of length 2L and width 2w considerably 
exceeding its thickness d, 2L > 2w 3> d, one may ap- 
ply first the field H az perpendicular to the plane of the 
strip, and then an oscillating in-plane magnetic field H ax 
across the width of the strip. This "shaking" leads to a 
homogeneous distribution of the perpendicular field H az 
over the sampled After this shaking process one keeps 
H ax = constant and applies the field H ay along the axis 
of the strip. 



Jc±a 




where D 2 = 1 — cos 2 ^ — -0) sin 2 9 w 1— sin 2 <^sin 2 9 (since 
either ^ a tt/2 or sin 2 9 « 1), cos 2 9 a H 2 z /(H 2 Z +H* y ), 
and the length a is determined by the sheet current J and 
cos#, 



J 2a d 

— cos 9 — — In — 
J r d \2a 




2a . d . . 
— arcsmh — . (46) 
d 2a 



We shall denote the solution of Eq. (|4"6")) as 2a /d — 
g(J cos9 / J c ). The function g(u) defined by u = 
garcsinh(l/(7) increases monotonically with its argument 
u, Fig. 6. Hence with increasing H ay , i.e., with decreasing 
cos 9, the length a decreases. When J is close to J c and 
cos 9 « 1, the length a tends to oo, while for J <C J c one 



has 2a jd ~ J cos9/J c C 1. A very good approximation 
valid at all u is^- 

. I (i + ,g 2 ) 1/2 + i 

g(u) w u/ In y , 

/ On 



9o(u) 



u + u l 



[24(1 - u)]i/2 



(47) 



The field of full flux penetration can be estimated from 

hy(0)=0, 



H° v = jc± a(Jl + ^-l), 



(48) 



where a is determined by Eq. (|46|) . When i? Qy <C Haz, 
one has cos 9 ~ 1, and the length a is almost independent 
of H ay . Thus, at such H ay the profiles <p{z), h x (z), and 
hy(z) — h ay practically do not change with increasing 
H ay - Fig. 5 shows that this property of the profiles, in 
fact, holds in the region H ay < H az /2 when J/J c = 0.5. 
However, if J is close to J c , the length a sharply depends 
on cos#, Fig. 6, and the width of this region shrinks. 
Solution (|45|) can be obtained as follows: We put 



(en ±x )" = 0, (en ±y )" = 0, 



(49) 



since it may be verified that the term proportional to 
ip' sin 2 9 in Eq. and the left hand side of Eq. fl30|) are 
small in the parameter J c / H az and hence may be omitted 
in the first approximation. Equations (|49p mean that 
enj_ x and en± y are linear functions of z, and thus each 
of them generally depends on two constants. However, 
taking into account boundary conditions (|43[) and the 
symmetry of the problem, one finds that the functions 
en± x and en± y are expressed by only one constant that 
coincides with enj_ y . If we denote this constant as a cos 9 
and use 



n± x 



cos ip 



sin <p cos 
D 



(50) 




FIG. 6: The function g(u) denned by u — garcsinh(l/<7) (cir- 
cles). Also shown is approximation (|47[l (solid line). 




d 
slab 



FIG. 7: Schematic picture of the strip. Shown is also a "slab" 
cut-out near the point x, see Sec. Ill C. 



we arrive at formulas (|45p . The profiles h x (z) and h y (z) 
follow from Eqs. (|26|) and (f27|) . and the constant a can 
be found from the condition 



J=Jy = 



d ^ 2 dz j c ± sin (p 

-d/2 



D 



(51) 



which is just Eq. (j46| . 

It is also instructive to write the electric-field com- 
ponents E x = erix x and E y = enj_ y explicitly. Using 
Eqs. (@SJ) and ([50]). we find 



E x = —noH ay z 
E y = [i H ay acos9 



(52) 

fX ° H n ayd cos 8g(J cos 9/ J c ). (53) 



The field £J a results from the tilt of a vortex line along the 
y direction when H ay is applied to the slab. Note that 
I -d/2 E x (z)dz — since the upper (z > 0) and lower (z < 
0) parts of the vortex move in opposite directions when 
the tilt occurs. On the other hand, E x is independent of 
z. This component of the electric field is due to a drift 
of the vortex as a whole in the x direction when H ay is 
applied to the sampled 

The above formulas for the slab with a current en- 
able one to reproduce a number of results for the 
vortex-shaking effect that were derived from geometrical 
considerations! 40 ' 41 In particular, the expression for (f(z) 
in Eqs. (f45|) . formula (j46|) . and Eq. (|53l) . in fact, coincide 
with Eqs. (4), (6) and (28) from Ref. l4ll in which the so- 
called longitudinal vortex-shaking effect in a thin strip 
was considered. To obtain the formulas for the vortex- 
shaking effect in a rectangular platelet^ - one should con- 
sider the slab with H az H ay , J c and with the total 
current J flowing at an arbitrary angle to the y axis, i.e., 
when J = (J x , J y ). The appropriate solution of the crit- 
ical state equations is still obtained from Eqs. (|4"9"]l . but 
now there is no more symmetry restriction on the z de- 
pendences of all the functions, and en± x and en^ y are 
expressed via two constants. Similarly to Eq. (fSTj) . these 
constants can be expressed via J x and J y , and the solu- 
tion thus obtained reproduces the appropriate results of 
Ref. IM 
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FIG. 8: The sheet current J(x) = \J y (x)\ in the strip to which 
first a large perpendicular magnetic field H az = 20 is applied 
and then an increasing longitudinal field H ay . The aspect 
ratio of the strip is 2w/d = 20. The magnetic fields are in 
units of J c — j c ±d. 



C. Third example: strip 

We now consider the third example of the general T- 
critical state. Let a thin strip fill the space |x| < w, 
\y\ < oo, \z\ < d/2 (d <C w), and be in a constant and 
uniform external magnetic field H az directed along the 
z axis, i.e., perpendicularly to the strip plane. The crit- 
ical current density j c ± is still assumed to be constant, 
and let H az considerably exceed J c = j c ±d so that at 
the initial moment of time, t = 0, the strip is in the 
fully penetrated Bean critical state. In other words, the 
magnetic- field profile H z (x) in the strip is described by 
the well-known function^&£ and one has J v (x) — J c for 
— w < x < and J y (x) — — J c for w > x > 0, where the 
sheet current J y is the current density integrated over the 
thickness d. At t > the magnetic field H ay is switched 
on in the y direction, and hence the applied field is tilted 
towards the axis of the strip. Note that the critical states 
in isotropic and anisotropic strips placed in inclined mag- 
netic fields were studied in Refs. [T3ll3Ii3IiMi3- How- 
ever, in all these papers the external magnetic field was 
tilted perpendicularly to the axis of the strip, the currents 
in the critical states were always perpendicular to local 
magnetic fields, and thus, the usual Bean critical states 
occurred in the strips. In the considered case the general 
T-critical states develop in the strip, and these states 
differ from the states of the second example in that the 
magnetic field H z and the currents J are not uniform in 
the x direction any more. 

Strictly speaking, the description of the magnetic-field 
tilt towards the axis of the strip reduces to solving a 
two-dimensional general T-critical state problem. How- 
ever, the smallness of the parameter d/w enables us to 
simplify this problem by application of the approach of 
Ref. 133 . Within this approach, we split the problem into 




x / w 

FIG. 9: The perpendicular magnetic field H z (x) caused by 
the sheet current of Fig. [8] The magnetic fields are in units 
of Jc 

two simpler problems: A one-dimensional problem across 
the thickness of the sample, and a problem for the in- 
finitely thin strip. Namely, we first interpret a small sec- 
tion of the strip around an arbitrary point x (see Fig. 7) 
as an "infinite" slab of thickness d placed in a perpendic- 
ular dc magnetic field H z (x) and in a parallel field H ay 
and carrying a sheet current J y (x). This is just the prob- 
lem that has been solved in Sec. Ill B. We then use the 
resulting electric field E y obtained for the slab, Eq. I]53p. 
as the local electric field E y (x) for an infinitely thin strip, 
to calculate the temporal evolution of the sheet current 
J y (x) and of the magnetic field H z (x) in this strip by the 
method of Ref. [631164 . The resulting equation for J y (x, t) 
can be written in the formi^ 

dJ y (x,t) _ _2_ f w _dv_ ( w 2 -u 2 ^ /2 dE y {J y ) 
dt ir no J- w u — x \w 2 ~x 2 ) du ' 

(54) 

where E y (J y ) is given by Eq. (|53j) . On determin- 
ing J y (x,t), the magnetic- field profiles are found from 
the Biot-Savart law. Since E y (J y ) oc H ay , we see 
again that the temporal dependence of the current 
and magnetic-field profiles is only a parameterization 
of their dependence on H ay , Sec. II A. It also fol- 
lows from Eqs. (|53[) and (|54f that these profiles de- 
pend on the parameters H ay , H az , d, w via the fol- 
lowing combinations: J y = J y (x/w, H ay / H az , P), H z — 
H z (x/w, H ay /H azi P) where we have introduced the no- 
tation P = (d/2w)H az / J c . Note that the considered 
critical state problem is similar to the problem of the 
longitudinal vortex-shaking effect in a thin strip4i The 
difference between the problems is that the magnetic field 
H ay now increases monotonically rather than oscillates 
about H ay — 0, and here we present results up to large 
values of H ay even as compared with H az . 

In Figs. 8 and 9 we show the profiles J(x,H ay ) = 
| J y (x, H ay )\ and H z (x, H ay ) that develop in the strip 
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during increase of the longitudinal field component H ay , 
i.e., when the applied field is tilted away from the z 
axis towards the strip axis y. The profiles J(x,H ay ) 
take the shape similar to the shape of the profiles in 
the longitudinal vortex-shaking effect^, and their mag- 
nitude decreases with increasing H ay . However, in con- 
trast to the shaking effect, this magnitude does not de- 
crease down to zero but tends to a finite limit that 
depends on the only parameter P = [d/2w){H az J J c ). 
Thus, at H ay 3> H az the current profiles J(x, H ay ) 
and the magnetic-field profiles H z (x, H ay ) reach nonzero 
limiting distributions. The existence of such limiting 
J(x) and H z {x) can be understood from the following 
considerations: At small cos#, if one neglects logarith- 
mic corrections, the electric field E y , Eq. (j53"j) . is pro- 
portional to H ay (d/2w)(J/ J c ) cos 2 8, and equation (j54")) 
has a solution with separable variables: J y (x, H ay ) = 
J c f(x/w)F(H ay ) where f(x/w) and F(H ay ) are some 
functions and cos 2 — H% Z /(H% Z + Hl y ). Inserting this 
form of J y (x,H ay ) into Eq. (|54p . we find that 



In F(H a 



-Parctan — — 



const., 



(55) 



i.e., at H ay — > oo the function F(H ay ) does not tend to 
zero. In other words, with increasing H ay the decay rate 
of J decreases so quickly that J does not reach zero even 



in the limit H ay — > oo. 



In Fig. 10 we show the magnetic moment per unit 
length of the strip, M z , in the general T-critical states 
developed during increase of H ay , 



M z {H ay ) = -2 x\J y (x, H ay )\ dx . 
Jo 



(56) 



Here J y (x, H ay ) is the solution of Eq. (|5"4")) . see Fig. 8. At 
the initial moment of time, when the usual Bean critical 
state occurs, we have M z (0) = —J c w 2 . The application 
of H ay leads to the relaxation of M z towards a saturation 
value M z (oo). The normalized moment M z (H ay )/M z (0) 
depends only on H ay /H az and the parameter P. All the 
curves of Fig. 10 may be well fitted by stretched expo- 
nentials of the form s + (1 — s) cxp[—p(H ay /H az ) q }, with 
some s and p depending on P and with q close to 2/3. 65 
However, an expression useful for all P is suggested by 
Fig. 10 (bottom), namely, 



M z {H ay 
M z (0) 



exp{-1.67 [Pavctan{H ay /H az )} ob5 } . (57) 



This expression, depicted in Fig. 10 (top and bottom) as 
dots, gives excellent fits to the numerical results for not 
too small P > 0.5, but even for smaller P it is qualita- 
tively correct and only slightly underestimates the exact 
M z at large H ay /H az . 

The saturation values s — M z (oo)/M z (0) are deter- 
mined by the above-mentioned limiting current profiles 
and depend on the only parameter P — {d/2w)H az / J c . 
These values obtained numerically are plotted in the in- 
set of Fig. 10 as circles, while the solid line in this inset 



~ 0.5 




1 H 1 ' 5 /H 2 

ay az 





M z H / M z (0) 




saturation 




\ value 



0.001 0.01 



Pxarctan(H /H ) 

v ay az' 



FIG. 10: The perpendicular component M z of the magnetic 
moment of the strip shown in Fig. 7, plotted versus Hay / H a z 
(top) and versus P ■ arctan(H ay / H az ) (bottom) for different 
values of the parameter P = (d/2w)H az / J c = 0.05, 0.1, 0.2, 
0.5, 1, 5. The dots show approximation (|57[) . The inset shows 
the saturation values s — M z (oo)/M z (0) calculated numeri- 
cally (circles) and their fit by Eq. (f58|) (solid line). 



is the following analytic approximation, 



s(P) w (0.5 - 0.5tanh{0.41[ln(P) - 0.5]}) ' 



(58) 



When the magnetic field H ay is switched on, not 
only does the z component M z of the magnetic moment 
change but there appears also a magnetic moment M y 
along the axis of the strip. This moment (per unit length 
of the strip) is defined by the expression: 



i-d/2 



/W rO 
dx / dzj x z, 
-w J-d/2 



(59) 



where the x component of the current density, j x = 
j c ± cos (p/D, can be found using solution (I45[) . With 
Eq. (|27p . formula (]59[) can be rewritten in the form: 

/w rd/2 
dx / dz[H ay -hy(x,z)]. (60) 
-w J-d/2 
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(d / 2w) H /J 

ay c 



FIG. 11: The longitudinal component M y of the magnetic 
moment of the strip shown in Fig. 7, plotted versus H ay /H az 
(top) and versus P ■ H ay /H az — (d/2w)H ay / J c (bottom) for 
the same P as in Fig. 10. Here M y (oo) = M° = -j c± d 2 w/2. 
Note that the curves for different P > 0.5 practically coincide 
if the magnetic moment is plotted versus P ■ H ay /H az . The 
dots indicate the fit (|62p for this universal curve. 



In other words, M y is the "expelled" flux in the y direc- 
tion. Inserting Eqs. (|45[) into this formula, we obtain 



My (H ay 

MO 



\ — h/T 

Jo w 



+ g 2 (u)-ug(u)], (61) 



where M® = —j c ±d 2 w/2 is the magnetic moment in the 
fully penetrated Bean critical state which occurs if the 
field H ay alone is applied to the strip, u = J cosO/ J c , 
and J = J(x, H ay ) is the current profile obtained from 
Eq. IJHJ, see Fig. 8. 

Figure 11 (top) shows the normalized magnetic mo- 
ment M y (H ay ) / M ay (oo) plotted versus H ay /H az for the 
same values of the parameter P as in Fig. 10. The 
saturation value M ay (oo) always coincides with M y = 
—j c j_d 2 w/2. Beside this, we find numerically the fol- 
lowing interesting result: If P is not too small, P > 
0.5, the normalized magnetic moment plotted versus 
P ■ H ay /Haz = (d/2w)H ay / J c is well described by the 
unique curve, Fig. 11 (bottom), 



At smaller P values fits of the form (62) are still possible, 
but with different fitting parameters. 

The results of this section describe the relaxation of 
J(x,H ay ) and H z (x,H ay ) to the limiting profiles and 
of M z (H ay ) and M y (H ay ) to their saturation values. 
According to Figs. 8-11 and Eqs. ([57]). (j6"2")l . this re- 
laxation mainly finishes at some H ay proportional to 
mm(H az ,2wJ c /d) [note that H az J c = dj c ± for 
Eq. (|5"3"|) to be valid and for the full flux penetration to 
occur in the initial state]. All these results can be veri 
fied in experiments similar, e.g., to those of Refs. 
However, we emphasize that in contrast to Refs. 



66 
66 



the magnetic-field component H az has to be switched 
on before H ay . This guarantees absence of flux-line cut- 
ting for not too large H ay , see Eq. (TM1) . If similarly to 



experiment 



,66,67 



the in-plane magnetic field is switched 



My (H ay ) 

My (oo) 



1 - exp{-4.8 [(d/2w)(Hay/J c )] 0Ai } . (62) 



on before H az , completely different critical states will de- 
velop. 



IV. CONCLUSIONS 

In this paper we point out how to calculate the general 
T-critical (cutting-free) states in an arbitrarily-shaped 
type-II superconductor when the applied magnetic field 
H Q slowly changes in its magnitude and direction. In 
accordance with the definition of the general T-critical 
state, it is assumed here that the external magnetic field 
changes in such a manner that flux cutting does not oc- 
cur in the sample. Our approach enables one to take 
into account the anisotropy of flux-line pinning and the 
dependence of the critical current density perpendicular 
to a local magnetic field, j c ±, on the longitudinal com- 
ponent of the current density jn . We also show that the 
variational principle recently propose d 44 i 45 > 46 cannot give 
the correct description of the general T-critical states for 
many situations. 

We analyze three examples of the general T-critical 
states, at least two of which may be investigated exper- 
imentally. In particular, we study a seemingly simple 
problem that has not been solved as yet, viz., we consider 
the critical states in a slab placed in a uniform perpendic- 
ular magnetic field H az and then two components of the 
in-plane magnetic field, H ax and H ay , are applied suc- 
cessively, Sec. Ill A. We obtain that one of the in-plane 
components of the magnetic moment, M x , becomes pos- 
itive with increasing H ay for any sign of M x in the initial 
state (i.e., at H ay = 0). This paramagnetic effect is due 
to the fact that the currents in the critical states are not 
perpendicular to the local magnetic fields. This effect is 
especially evident when H az is of the order of the self- 
fields of the slab. 

In the other example, we analyze the general T-critical 
states in a long thin strip placed in a perpendicular mag- 
netic field H az which then tilts towards the axis of the 
strip y, Sec. Ill C. When H ay , the axial component of 
the applied magnetic field, increases, the magnetic-field 
and current profiles across the width of the strip tend to 
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limiting profiles, and the components of the magnetic mo- 
ment, M z and M yi reach saturation values. The limiting 
profiles and the saturation value M z (oo) for M z (H ay ) are 
determined by the only parameter P = (d/2w)H az / J c 
where d and 2w are the thickness and the width of the 
strip, respectively, and J c — dj c ±. If P is not too large, 
P < 5, the limiting profiles and M z (oo) noticeably differ 
from zero, while at P > 5 they become very small and 



practically vanish. The saturation value for M y is always 
equal to M° = -j c ^d 2 w/2. 
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